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The aim of this paper is to give a precise estimate on the tail prob- 
ability of the visibility function in a germ-grain model: this function 
is defined as the length of the longest ray starting at the origin that 
does not intersect an obstacle in a Boolean model. We proceed in two 
or more dimensions using coverage techniques. Moreover, convergence 
results involving a type I extreme value distribution are shown in the 
two particular cases of small obstacles or a large obstacle-free region. 

1. Presentation of the model and results. In [18] G. Polya intro- 
duced the question of the visibility in a forest in a discrete lattice case as 
well as in a random case. He first treated the problem of a person standing 
at the origin of the regular square lattice of M? , when identical trees (discs 
with constant radius R) are situated at the other points of the lattice. In this 
framework he showed that in order to see at a distance r the radius R should 
be (asymptotically when r is large) taken as 1/r. More recently V. Jankovic 
gave in [11] an elegant proof of a detailed version of this result. The random 
case studied by G. Polya was the one of the visibility in one direction: we are 
here interested in the global solution to this problem considering all direc- 
tions simultaneously. The spherical contact distribution which can be seen 
as the infimum of the visibility over all directions has been intensively used 
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for a geometric description of random media (see e.g. [1, 5, 10, 14, 16, 19]). 
In comparison, the total visibility, i.e. the supremum of the same function, 
has been rarely studied in the litterature. The work of reference is due to S. 
Zacks and is strongly motivated by military applications ([25], see also his 
work with M. Yadin [24]). However his interest was mainly focused on the 
probability that given points could be seen and not on the total visibility. 
Very recently, the visibility problem has been investigated in the hyperbolic 
disc by I. Benjamini et alt. [2]. In particular, the authors show the existence 
of a critical intensity for the almost sure visibility at infinity. In this connec- 
tion, one of the consequences of our work will be that with probability one, 
we can see only at a finite distance in the Euclidean space M (see Proposi- 
tion 4.1 and also Remark 4.1 concerning the possibility to see at infinity). 

In this paper, one of our goals is to present new distributional proper- 
ties of the total visibility in order to develop a future use of this indicator 
for the study of porous media and more particularly in forestry. Potential 
applications concern the optimization of directional logging of trees, the 
measurement of competition level between growing trees in forest dynamics 
or even an estimation of the light transmission through the canopy of a tree. 
In such context, the total visibility seems to have an important role to play 
even though it is understood that in some particular cases, another quantity 
of interest could be the mean of the visibility in all directions. 

The model is the following: consider a Boolean model (see [15, 23]) with 
random almost surely diameter-bounded convex grain K with law [i based 
on a Poisson Point Process X with intensity measure the Lebesgue measure 
on R , d > 2. Define £> the occupied phase of this model, 

O = (J (x®K x ), 

where (K x ) x ^x are independent identically distributed copies of K, indepen- 
dent of X. We condition this model by the event O D where O is the origin 
of M. d . In particular, it has a positive probability equal to exp(— E[Leb^(K)]). 
We then define the visibility in the following way: 
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Fig 1. Directional visibility 

Definition 1.1. Let u be a unit vector in M. d , the visibility in direction 
u is defined as 

V(u) = inf{r > : ru G O}, 
the total visibility is defined as 

23 = sup V(u). 

||u||=l 

For a convex body K and u E S d_1 , we define the width of K in the 
direction v orthogonal to u as 

W U {K) = sup {y - x, v). 

x,y£K 

The mean width of K is denoted by W(if) and satisfies 
W(K) = ^-| ^(^d^H, 

where is the uniform measure on S rf_1 and a;^ is the d-dimensional 
Lebesgue measure of the unit-ball in M. d . 

It is well known that the law of the directional visibility is exponential, 
indeed 

Lemma 1.1. For each unit direction u in M. d one has 
P(V(u) > r) = exp(-rE[W(if)]). 

The aim of this paper is to give precise estimates on the tail probabil- 
ity for the visibility function in all directions: in section 2 we present the 
general method of coverage processes used throughout this paper. In sec- 
tion 3 we give sharp upper and lower bounds in the same exponential order 
exp(— Const. r) in dimension two and in the two cases of circular obstacles 
and of more general rotation-invariant random obstacles. In higher dimen- 
sions results on coverage processes are more sparse, thus the bounds pre- 
sented in section 4 for dimension d > 3 are rougher. Section 5 is devoted 
to two similar convergence results for the asymptotics of the visibility with 
small obstacles, and when the spherical contact length is conditionned to be 
large. Both results state a convergence in law towards a Gumbel distribu- 
tion, are valid for any dimension d > 2 and are based on an extension of a 
result of [12]. 
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2. Random coverage of the circle and the sphere and visibility. 

The visibility up to length r may be blocked only by those obstacles that 
intersect Bd(0,r), the number N r of those obstacles is Poisson distributed 
with parameter 

g d (r, K) = K[Leh d ((B d (0, r)®K)\ K)}, 

where K = {—x : x G K}, and those obstacles Ki, i G \1, . . . ,N r }, are 
independent and identically distributed with the same law. Each of them 
projects a shadow of solid angle Si on the sphere r§ rf_1 = dB d (0,r), which 
is a random cap Cf. 

(2.1) d = {(r, u) G M+ x S d_1 : 3s < r s.t. (r, u) G Ki}, 
and the solid angle Si is equal to 

(2.2) Si = {ueS d ~ 1 : (r,u)eC t }. 
We have thus the following ansatz: 

Ansatz 2.1. The visibility QJ is greater than r if and only if the solid 
angles Si, i G {1, . . . , N r } do not cover the sphere 

This equivalence links coverage properties of the sphere with our initial 
problem: this problem of random coverings has been quite intensively stud- 
ied in the literature in the two-dimensional case, see for instance [6, 7, 13] 
in the context of Dvoretzky covering, and [20, 21, 22, 24] for a more general 
approach. The properties of coverings in higher dimensions are less known, 
let us cite the works of S. Janson ([12] and other papers) dealing with some 
asymptotic properties for coverage processes with small caps. 



Let us remark that each obstacle Ki is distributed according to the law 

g d (r, K) _1 l (:re ^ )nBd(Oir) - ! .0 dxdn(K). 

We denote by v T the distribution of the associated solid angle Si . The prob- 
ability measure i> r is naturally invariant under rotations. We call P(p r ,n) 
the probability that n independent solid angles distributed as v r cover the 
sphere. With Ansatz 2.1 the following result becomes straightforward: 

Proposition 2.1. For every r > 0, 

(2.3) P (53 > r) = exp (-g d (r, K)) £ (9d(r,K)) n _ p 
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When in dimension two, for every probability measure v on [0, 1], P(v, n) 
denotes the covering probability of the circle with perimeter one by n i.i.d. 
isotropic arcs of ^-distributed length. We shall denote by v r the probability 
law of the random arc shadowed by an obstacle intersecting the disc -E?2(0, r). 
When K = E>2(0,R) for a constant radius R, one has for instance the fol- 
lowing values for the characteristic data of proposition 2.1: 

(2.4) g 2 {r,K) = n(2rR + r 2 ), 

df (u) = ^T^ 1 [0-^ct aa (f)]W >< 

. ,„ . sm(iru)(R 2 + r 2 cos(27ra)) \ 
rsm{2TTU) H 

J R 2 — r 2 sin 2 (-7rn) / 
ttR 2 cos(7rn) 

(2 - 5) + 7r~T^ 8 -^U] {u) sTn^) • 

Remark 2.1 (The two-dimensional case). Let us recall that the covering 
probability P (u, n) has been computed in dimesion two by Stevens in the case 
of arcs with deterministic lengths [22] and Siegel & Hoist in the general case 
[21 J: for every probability measure v on [0, 1] and n £ N* ; 
(2.6) 



P{v,n) 




k 



=i 



i=i j ° 



-k 



dA fc (n) 



where F u is the cumulative distribution function of v and k G N*, is 
the normalised uniform measure on the simplex {(x\, ■ ■ ■ ,Xk) G [0, l] fc ;xi + 
\-x k = l}. 

We can consequently substitute P(v,n) by its expression (2.6) to get an 
explicit formula for P (5J > r), r > 0. Nevertheless, it seems more or less 
intractable for doing some asymptotic estimations. 

3. Sharp asymptotics in dimension two. In dimension two it be- 
comes possible to give sharp estimates for the tail probability, for instance 
we have from lemma 1.1 the following lower bound: 

(3.1) P (53 > r) > P (V u > r) = exp(-rE[W(iT)]), 

in this section we shall give two sharper lower bounds and an upper bound: 

• When obstacles are fixed discs with constant radius R, for each /j, G 
(0, 2/R) there exists a function e M converging to as r — > oo such that: 

(3.2) P(5J > r) > firexp(-2Rr) (1 + e M (r)) , 



imsart-aap ver. 2009/02/27 file: submission-visibility.tex date: May 29, 2009 



G 



P. CALKA ET AL. 



for r large enough; 
• When obstacles are discs with random radius R, we have for r large 
enough 

P(53>r) < 2(ir(2rE[R] + r 2 ) + 2) x 

(3.3) exp (-ir(2rE[R] + r 2 )m r ) , 

P(9J > r) > 27rm r (2rE[R] + r 2 ) exp (-Trm r (2rE[R] + r 2 )) + 

(3.4) exp (-27rm r (2rE[ J R] + r 2 )) , 

where m r is the mean of law v r , satisfying 

, . 2ELR1 

(3.5 m r 

7rr 

From (3.3), (3.4) and (3.5) we obtain directly the following theorem: 

Theorem 3.1. For the Boolean model in dimension two with random 
discs, the asymptotics of the visibility is given by 

lim - log P (53 > r) = -2ELR1. 

r— »+oo r 

The estimate (3.2) is obtained without the use of the previous section, by 
more direct considerations based on the estimation of the visibility in a finite 
number of directions. Details of the proof are postponed to the appendix 1. 

In the following lines, we shall prove estimates (3.3) and (3.4), as well 
as a generalisation of theorem 3.1 for generall convex shapes. Most of the 
arguments rely on comparison results for covering probabilities of the circle. 



3.1. Lower bound via comparison of coverage probabilities, estimate (3.3). 
In this subsection and the next we shall follow the main steps of a previous 
work [4] that dealt with the circumscribed radius of the typical Poisson- 
Voronoi cell: a comparison result on covering probabilities states that if two 
probability measures (the laws of the arcs) are comparable in some sense, 
then the coverage probabilities are also comparable. Up to now the ordering 
induced by the concentration around the mean has been the main criterion 
for comparing covering probabilities [4, 20] but the convex ordering (which 
is implied by the previous ordering) is in fact enough to deduce the required 
inequalities. 

For the special case of random discs, the computation of the parameters of 
equation (2.3) is done in the following way. Let us denote by \x the law of the 
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radius of the random discs of the Boolean model. The law of the normalised 
lengths of the shadowed arcs on the circle is denoted by u r , this law is the 
image of the couple (R, U) with law /j, <g> U(0, 1) by the map (R, U) t— > I: 



(3.6) t 



(I 1 r 

— arcsin — - if < U < ; 

y/l + I (2 + £) U " r + 2ii 

I arccos £ + ( 2 + £) U _ if <U<1. 
* 2Jl + ^(2 + ^) [7 r + 2i? 



As r tends to infinity, we obtain £ ~ (1/tt) arcsin(i?/(r\/?7)) so that the 
asymptotics of the expectation becomes m r := ~E[l] ~ 2E[i?]/(-7rr). 

We shall use here the convex domination of measures, let us first remark 
the following: for any probability measure v on [0,1/2] with mean m and 
any convex function / : [0, 1/2] — > R, we have 

(3.7) J f dv < (1 - 2m)/(0) + 2m/(l/2), 

which means that v < cv (1 — 2m)5o + {2m)8\/2 (where < cv denotes the 
usual convex order [17]). 

It is a consequence of the proof of theorem 13 in [4] that the convex order 
implies the order of the covering probabilities: if \i\ and \ii are two proba- 
bility measures on [0,1/2] such that fj,% < cv fj-2 then P(^i,n) < P(/J,2,n). 
Inserting the inequality P(v r ,n) < P((l — 2m r ) + 2m r 5i/2,n) in (2.3) for 
every n £ N, we get that 



' (9J > r) > exp (-ir(2rE[R] + 



r 2 ) I x 



(3 8) E ('i"^' ! )'- ( i _ P ( (1 _ 2mr)a „ + 2mr il/2 

n>0 



It follows from ([4], Corollary 1) that 

1 - P ((1 - 2m r )5 + 2m r 5 1/2 , nj = 2nm r (l - m r ) n ~ l + (1 - 2m r ) n . 

Inserting that result in (3.8), we obtain that 

P (53 > r) > 2vrm r (2rE[i?] + r 2 ) exp (-vrm r (2rE[ J R] + r 2 )) + 
(3.9) exp (-27rm r (2rE[i?] + r 2 )) . 
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3.2. Upper bound via comparison of coverage probabilities, estimate (3.4)- 
By Jensen's inequality, we have that v r > cv 5 mr . Consequently, the same 
argument as for the lower bound shows that P{y r ,n) > P(5 mr ,n) for every 
n S N. Inserting this inequality in (2.3), we have 
(3-10) 

P (!U > r) < exp (-n(2rE[R] + r 2 )) £ frfoE[fl] + r 2 ))" _ p 

For the estimation of (1 — P(5 a ,n)), Shepp obtained a basic inequality 
[20] which holds for a £ [0, 1/4] and n G N* and is easier to use than Steven's 
explicit formula: 

2d - a ) 2n 

(3.11) l-P(S a ,n)< 



j£(i_ _t)»dt+(! -a)(l-2o)»' 



A straightforward consequence of (3.11) is that for every a G [0,1/4] and 
n E N, we have 



(3.12) l-P(<y oj n)<2(n + l)(l-o 



In particular, for r sufficiently large, the mean m r is in the interval [0, 1/4] 
so the equality (2.3) combined with (3.12) leads us to 



IP (9J > r) < 2(vr(2rE[ J R] + r 2 ) + 2) exp(-vr(2rE[ J R] + r 2 )m r ). 
It remains to use the estimation on the mean m r to get that 

(3.13) limsup- logP(2J > r) < -2~E[R]. 

r— >+oo T 

(3.9) and (3.13) now complete the proof of Theorem 3.1 for random discs. 

3.3. The case of general convex shapes. Let K be a random convex body 
of M? containing the origin, which is supposed to be invariant under any 
rotation and is such that its diameter is bounded almost surely by a constant 
D > 0. For instance, /C can be the image of a deterministic convex body 
by a uniform random rotation. By the rotation-invariance of K, we have 
E[W U (K)]=E[W(K)\. 

Theorem 3.2. For the Boolean model with random rotation-invariant 
grains distributed as K , the asymptotics of the visibility is given by 

lim - log P (53 > r) = -E[W(K)]. 

r- ■> t oc r 
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Proof. The lower bound is obtained by lemma 1.1. It remains to show that 
limsup-logP(9J > r) < -E[W(K)]. 

r— >+oo f 

In order to do it, we need an intermediary geometric lemma whose proof is 
postponed to the appendix 2: 

Lemma 3.1. Let L be a convex body of diameter bounded by D containing 
the origin. We define $(ru + L) as the angle of vision of (ru + L) from O. 
Then 

lim r$(ru + L) = W U (L), 
r— >+oo 

the limit being uniform over all unit-vectors u and all such convex bodies L. 

As in the proof of Proposition 2.1, the event {23 > r}, r > 0, can be seen 
as the uncovering of the circle C(0,r) = rS 1 by the 'shadows' produced by 
the obstacles (x © K x ) such that (x © K x ) n B 2 (0, r) / 0. 

Let e > 0. By Lemma 3.1, let us fix r e > such that for every x such 
that ||x[| > r e and every convex body L (with a diameter bounded by D), 
we have 

(3.14) ||x||*(x + L) > (W x/M (L)-e). 

Then for r > r e + D, the probability of uncovering the circle C(0, r) is greater 
if we only keep the shadows produced by the obstacles (x © K x ) such that 
r e < \\x\\ < r — D. In that case, such a shadow is a random rotation-invariant 
arc on the circle C(0, r) whose normalised length is (27r)~ 1 ^I , (xffii ; C a ;). Let us 
denote by r\ r the mean of (2tt)~ 1i ^(Z Q) Kz) when Z is uniformly distributed 
in B 2 (r — D)\ £?2(r e ) and Kz is independent from Z and distributed as IT. 
Following the method already used to obtain the upper-bound (3.10), we 
have 

P(© > r) < e-^- D ^ £ *((r-g) a -rg) (1 _ p{§ )y 

n=0 

Forecasting that r\ r will be small enough, we may apply the inequality (3.12) 
in order to obtain that 

(3.15) F (2J > r) < 2(vr((r - D) 2 - r 2 £ ) + 2) exp(-vr((r - D) 2 - r 2 )rj r ). 
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Let us now estimate the mean rj r : using (3.14), we get 



2tt 



Ek 



> 



2tt rr — D 

/ ^(pu e ®K)p 

Jr e 

2 \w e (K)-e)&9 



dp 



vr((r - D) 2 - r 2 ) 
r - D dp 



de 



7r((r - D) 2 - r 2 ) 



r— *+oo 



(E[W(tf)]-e) 
E[W(iT)] - e 



r-D 



dp 



e vr((r - D)2 - r|) 



71V 



Inserting this last result in (3.15), we obtain that 

limsup-logP(Q3 > r) < -(E[W(iT)] -e). 

r— >+oo r 

When e goes to 0, we get the required result. 



□ 



4. Rough estimates in dimension greater than three. The prob- 
lem of maximal visibility in a Boolean model is investigated in M. d with 
deterministic radii R x = R, x £ X. The obstacles are balls of determinis- 
tic radius R. The same connection between the distribution of 03 and the 
non-covering of the sphere by random circular caps occurs. What prevents 
us from obtaining the analogue of Theorem 3.2 is that the calculation of the 
probability to cover the sphere with caps of random radii is not known. We 
have to restrict ourselves to coverings of the sphere with caps of determin- 
istic radii. This explains that the following result is weaker than Theorem 
3.2: 



Proposition 4.1. In dimension d > 3, we have 

liminf- log P (53 > r) > -uj^iR 4 - 1 

r— >+oo f 

and 

limsup- log P (53 > r) < --u d ^ 1 R d ~ 1 . 

r— »+oo f d 

Proof. As in the two-dimensional case, the lower-bound is obtained by 
considering the visibility in a fixed direction. 

Let us focus on the upper-bound: the maximal visibility is larger than 
r > if and only if the shadows produced by the obstacles on the sphere 
centred at the origin and of radius r do not cover that sphere. The concerned 
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balls are those such that their centres are at distance p £ [R, R + r] from the 
origin. Since we look for an upper-bound of a probability of non-covering of 
the sphere by random circular caps, we can take less and smaller caps. For 
sake of simplicity, we only keep the shadows produced by the balls with a 
centre at distance p G [R, r] . 

For such a ball, it comes from (3.6) that the angular radius of its shadow 
on the sphere is at least arcsin f^rjj which is bigger than —• 

In conclusion, the probability that the maximal visibility is greater than 
r is lesser than the probability of non-covering of the unit-sphere by a Pois- 
sonian number (of mean tUd{r d — R d )) of circular caps of angular radius ^. 
Upper-bounds for covering probabilities of the unit-sphere have been pro- 
vided by Gilbert [8] in dimension three, Hall [9] in any dimension when the 
unit-sphere is replaced by the unit-cube and more recently by Biirgisser, 
Cucker & Lotz [3]. A very minor consequence of Theorem 1.1 of this last 
work is the following: let P(f, n) denote the probability to cover the unit- 
sphere with n random circular caps which are independent, with uniformly- 
distributed centres and with a fractional area of /. Then 

(4.1) lim -log(l-P(/,n)) = log(l-/). 

n— >+oo n 

In particular, the fractional area occupied by a circular cap of angular radius 
R/r is 



duJd Jo r-^+oo doj^r 

Consequently, we have 

P (QJ > r) < e"^-^) £ Mr d -R d )) n {l _ p {f ^ n)) 

n=0 n ' 

and a direct application of (4.1) provides that 

logP(S3>r) < -LO d r d f r ~ -^±± R d-\. 

rw+oo r^+oo d 

This completes the proof of Proposition 4.1. □ 

Remark 4.1. Proposition 4-1 implies that the total visibility 53 is finite 
almost surely. Nevertheless, when the intensity measure of the underlying 
Poisson point process is of the form (in spherical coordinates) r a ~ l drdad{u), 
a G R, it can be shown in the same way that the visibility at infinity exists 
with positive probability as soon as a < 1. 
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5. Small or distant obstacles: convergence towards the law of 
extreme values. When the size of the covering objects becomes smaller 
and the number of objects grows at the same time accordingly, S. Janson [12] 
showed in very general setting that a particular scaling yields a convergence 
towards the Gumbel law. We shall use this kind of result in two contexts 
below: 

• the asymptotics of the visibility when the obstacles are small (or equiv- 
alently when the intensity of the centres is small) ; 

• the study of the visibility when there exists a large region around the 
origin with no obstacle at all. 

5.1. Small obstacles. In this subsection, the radius R of the obstacles 
will be deterministic but no longer constant. For sake of clarity, we will 
denote by the visibility when the obstacles are discs of radius R > 0. 
We aim at giving the asymptotic behaviour of the visibility when the size of 
the obstacles goes to 0. Let us define the quantity 

(5.1) U = u d ^R d - x V3 R + d{d - 1) \og(R) - 2{d - 1) log(- log(i?)) - K d 
where 



f d 2(d-l)( d _ 1) 3(d-l)-l T U_l\ 2d 2 

^ (d-l)!vr^4^ i 2 M -3r( 



V 2 

Theorem 5.1. When R goes to 0, the quantity £r (provided by (5.1)) 
converges in distribution to the extreme value distribution, i.e. for every 
u£R, 

limP(£fl<u) = exp(-e- u ). 
-/t — >{j 

Proof. The proof relies essentially on the application of a result due to Jan- 
son (Lemma 8.1. in [12]) about random coverings of a compact Riemannian 
manifold by small geodesic balls. 

As before, we exploit the connection between the cumulative distribution 
function of QJr and the probability of covering the sphere with circular caps: 

did- l)log(R) 2(d - 1) log(- log(.R)) 

V tod-i R ^>d-i R 

K d + u 1 



= P(the sphere of radius f(R) is covered 
(5.2) by circular caps coming from the obstacles) 
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where 

_ djd - 1) logCg) 2(d-l)log(-log(fl)) + 1 

Let us focus on this covering probability: the concerned obstacles are those 
such that their centres x are at distance p S (R,R + f(R)). Their number 
is Poisson distributed, of mean ojd {{R + f(R)) d - R d ). The set of all a;/||x||, 
where iGXfl [B(0, R + f(R)) \ B(0, R)], is a homogeneous Poisson point 
process on the unit-sphere of intensity ((f(R) + R) d — R d )/d. 

The induced shadow of each of theses obstacles is a geodesic ball on the 
unit-sphere of angular radius equal to arcsin(i?/p) if p € (R, \/R 2 + f(R) 2 ) 
and equal to arccos((f{R) 2 + p 2 -R 2 )/2f(R)p) for p e yR 2 + f{R) 2 ,R + r) 
(see (3.6)). In particular, it can be verified that the normalized geodesic 
radius Or of this circular cap satisfies that 

(5-3) — Or l[i,+oo)(«)-^T du - 

where clr = ^— >^ 0. Consequently, the required covering probability in 
(5.2) is the probability that the unit-sphere is covered by a Boolean model 
on the sphere of intensity Xr = iil&tlH ~ R ) suc h that the geodesic balls 
have i.i.d. radii distributed as Or (with Or satisfying the convergence (5.3)). 

It only remains to verify that all the hypotheses of Lemma 8.1. in [12] are 
satisfied (§ rf_1 being a {d — l)-dimensional Riemannian manifold): 

• the only notable difference is that we should not have ^@_R converg- 
ing in distribution but have it of fixed distribution for any R > 0. 
Nevertheless, the proof of Lemma 8.1. in [12] relies essentially on con- 
vergence results [(7.15), (7.20), ib.\ which also work in this context 
without any changes; 

• the moments of order {{d — 1) + e) of the limit distribution obtained 
in (5.3) are finite for every e 6 (0, 1). Moreover, the moment of order 
{d — 1) is d and the moment of order (d — 2) is d/2; 

• the constants b and a defined in ([12], Lemma 8.1) can be calculated: 

■ au 



dwd Ji u- LO d 



and 



a 




,x d-l 
d 



d d ~ 
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• The convergence (8.1) in [12] is satisfied: 




lim \ ba d R 1 dw d \ R + \og(ba d R 1 )+ 



(d-l)log(-log(6a|- 1 ))-log(a)} 



= u. 



Consequently, the proof of Theorem 5.1 is complete. 



□ 



Remark 5.1. The same type of method and result should also occur in 
dimension two when the discs are replaced by rotation-invariant i.i.d convex 
bodies. 

Remark 5.2. In any dimension, the result could be extended to radii of 
the form R x = eU x , x G X, where e goes to and the U x are i.i.d. bounded 
random variables. 

5.2. Conditioning by a large clearing. We define here S the clearing ra- 
dius as 



The distribution of S is called the spherical contact distribution. This section 
aims at estimating the distribution of the maximal visibility 53 condition- 
ally on S. In particular we show that when S is large, 53 is asymptotically 
equivalent to S (see Theorem 5.2) and we estimate precisely the difference 
(53 — S) via an extreme value result (see Theorem 5.3). 

A first estimation based on techniques similar to the proofs of (3.3) and 
(3.4) provides the following result: 

Theorem 5.2. For every a G (0, 1), we have 



Proof. Let us fix r > 0. Conditionally on {S > r}, the process of cou- 
ples (x,R x ) is a Poisson point process on R 2 x R + of intensity measure 

As in Proposition 2.1, 53 is greater than r + u, u > 0, if and only if the 
circle C(0,r + u) is not totally hidden by the 'shadows' of the obstacles. 
Moreover, the discs B2(x, R x ) which produce a non-empty shadow are those 
which satisfy r < ||sc|| — R x < (r + u). The formula for the length of the 
shadow depends on whether < \f[r + u) 2 + Rj, or not (see equalities 
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(3.6)). Consequently, if we only consider the shadows produced by the discs 
B<2(x,R x ) such that 



(r + R x < \\x\\ < J (r + u) 2 + R 2 ) and u > \/r 2 + 2rR x — r, 

then the probability of not covering the circle is greater. 
When u > \/r 2 + 2rR* — r, the number of such discs is Poissonian, of mean 
tt(u 2 + 2ru — 2rE(R)). Moreover, for these discs, the length of the shadow 
decreases with ||x|| and is minimal when ||x|| = \/Jr + u) 2 + R 2 , equal to 
L> min = (1/vr) &vcsm(R x / ((r + u) 2 + R x )) (see (3.6)). In the sequel, we denote 
by v r ,u the distribution of L min and Tn r u its mean. 

In conclusion, we have proved the following inequality: for every u > 
\J r 2 + 2rR* — r, 

F (9J > r + u\L > r) 

< e -.(. 2+2 _^(«)) g (^ + 2ru-2^(fl)))" _ 
n=0 

(5.4) < e-^+a™-*^) g ^ 2 + 2rn ; 2r ^^"(l - P(«5 mr>a! n)). 

n=0 

In particular, when u = r a , < a < 1, we have 

mrrQ ^ I^CH) an d vr(u 2 + 2ra- 2r J B(i?)) ~ 2vrr 1+Q 

' r— >+oo 7f 7" r— >+oo 

Using the inequality (3.12) and inserting the two previous estimates in (5.4), 
we obtain the required result, i. e. 



P (9J > r + r a \S > r) = O (V 2i?(fl)rQ 



Finally, it remains to study the distribution of 2J conditionally on S = r. We 
remark that conditionally on {S = r}, the process is the same as in the case 
of the conditioning on {S > r} with a supplementary random disc B%(Y, Ry) 
such that Ry is /i-distributed and conditionally on Ry, Y is uniformly dis- 
tributed on the circle C((D,r + Ry)- Since there is one more obstacle, the 
maximal visibility must be lesser than in the case of the conditioning on 
{S>r}. □ 
Theorem 5.2 implies that the difference (5J — S) is negligible in front of S 
but we can get a far more precise three-terms development in the following 
way: for every r > 0, we denote by 23 r a random variable distributed as 53 
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when the Boolean model is conditioned on {S > r}, i.e. on not having any 
grain at distance lesser than r from the origin. Let us define the quantity 

(5.5) Vv = Wd-i^iO (9J r -r)-(d-l) log(r) - (d- 1) log(log(r)) - K' d 

where 




d-2 

2\d-l 



d-l\ 



^d-iE^R ) 



doJd 



For every i £ R, we have 



(d - 1) log(r) - (d - 1) log(log(r)) - if < i|5 > r) . 

Theorem 5.3. When r goes to oo, the quantity ijj r (provided by (5.5)) 
converges in distribution to the extreme value distribution, i.e. for every 
t G R, 

lim P N) r <t) = exp (-e~ l ) . 

Proof. As previously in the case of small obstacles, the proof relies essentially 
on the application of a result due to Janson (Lemma 8.1. in [12]) about 
random coverings of a compact Riemannian manifold by small geodesic balls. 
Let us consider the quantity 



/(r) = ^X(R^) log(r) + «Je~(r^) log(log(r)) + 



such that (ip r < t) <^=^> (53 r — r < /(r)). The connection with a covering 
probability is the following: 

P (Vv < u ) = P(the sphere of radius (r + f(r)) is covered 

by circular caps coming from the obstacles). 

It remains to investigate asymptotics of this covering probability: the con- 
cerned obstacles are those such that their centers x are at distance p G 
(r+R, r+f(r)+R). Their number is Poisson distributed, of mean uv£^[(i?+ 
f{r)+r) d -(R + r) d }. 

The induced shadow of each of theses obstacles is a geodesic ball on the 
unit-sphere of angular radius equal to: 
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• aicsm(R/p) if p G (R + r, \/i? 2 + (/(r) + r) 2 ), 

• arccos(((/(r)+r) 2 + p 2 -i? 2 )/2(/(r)+r)p) if p > ^i? 2 + (/(r) + r) 2 
and p < R + r + /(r)) (see (3.6)). 

In particular, it can be verified that the normalized geodesic radius r of 
this circular cap satisfies that rQ r p. In the rest of the proof, we will use 
the quantity a r = 1/r in order to be as close as possible to the notations of 
Janson's lemma. 

Consequently, the required covering probability is the probability that the 
unit-sphere is covered by a Boolean model on the sphere of intensity 

\ r = ^[(R + f(r)+r) d -(R + r) d ] ~ r d ~ l f{r), 

such that the geodesic balls have i.i.d. radii distributed as r . 

As in the proof of Theorem 5.1, we verify that all the hypotheses of Lemma 
8.1. in [12] are satisfied: 

• the moments of order {{d— 1) +e) of the limit distribution p are finite 
for every e > 0. 

• the constants b and a defined in ([12], Lemma 8.1) can be calculated: 

b _ io d . 1 E^(R d - 1 ) 
du d 

and 

a ~(d-iy.[ r(|) ) ' e^r*-^ 

• The convergence (8.1) in [12] is satisfied: 

lim {bai^duJdXr + log(6a^ 1 )+ 

i — >+oo L 

(d - 1) log(- log(6a^ 1 )) - log(a)} = t. 

Consequently, the proof of Theorem 5.1 is complete. □ 

Appendix 1: proof of the lower bound via direct computation, 
estimate (3.2). It is quite reasonable to try to obtain directly a lower 
bound on the tail probability, let us explain the sketch of the proof: the 
visibility is greater than r if and only if there exists a direction in which one 
can see farther, so that if one discretises the circle 9^(0, r), one could argue 
that there exists one of those directions such that the visibility is greater 
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than r, the number of those directions is 2/wr, hence the order r exp(— 2rR). 
We shall make this statement more rigorous below. 

Let us take £ G (0,2/ R), we define N r as the integer part of (s and 
9 r = 27r/(£r), and define the points A k>r = (r, kO r ) for k € {0, . . . , N r — 1}. 
We see easily that if we define 

G k , r = (B 2 (0, R) © [0, A ktT ] ) \ B 2 (0, fl) , 

then for r > R and fei 7^ /c 2 one has 

Gfc!,r nG fc2 ,r C B 2 (0,r). 

The sets G k , r are called 'fingers', we denote by p r = Rj sin(# r /2) the max- 
imal norm of a point belonging to the intersection of two fingers, we shall 
denote by E^^ the intersection Gk t r H -62(0, p r ) and i^,. = Gfc ir \ E^ r . 

We will assume from now on that r is large enough. If at least one of those 
points Ak, r is not shadowed by the discs intersecting ^(0, r), the visibility 
53 is greater than r: hence the probability of this event is greater than the 
probability that one of the 'fingers' Gk, r in figure 2 does not contain a point 
of X. 




Fig 2. The points Ai )T and their associated 'fingers'. The point At, r ts visible if and only 
if no point of the process X belongs to Gk,r = Ek.r U Fk, r - 

We have 

Leb 2 (F kir ) = 2R(l - k)s + 0(1), 

where k = -RC7T -1 , so that 1 — k > 0, and 0(1) means a bounded function. 
Let us define the following events 

V k;r = {XnG kiT . = 0}, 

Zk, r = {x n Eh >r = 0} , 
w Kr = {xnF kir = 0}. 
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We have 

Vk,r = Zk,r H 

and we want to evaluate P (U/£o 1 Vk,r\ > this is equal thanks to Poincare's 

formula to 

(5.6) 

/N,— l \ N, — 1 N r 

P U V k,r )= E P + E(" 1 )^ 1 E P (*W n • • • n *W) • 

V fc=0 / k=0 1=2 ki<-<kt 

We shall prove that the dominating term in this expansion is the first one: 
it rewrites as 

N r -l N r -l 

F(V k>r ) = exp(-Leb 2 (Sfc,rUF feir )) 
k=0 k=0 

= N r exp(-2Rr). 
Let us consider £ > 2 and < k\ < ■ ■ ■ < kg < N r , we have: 

P (V kur n • • • n v ki>r ) = p ( f) z fci , r n f| w^,^ , 

\i=i i=i ) 

where the events Hi=i ^fc<,r! W klt r, ■ ■ ■ j Wfc^.r are independent as the sets 
Ui=i -^fci.r an d i^fc^r, • • • , i^, r are disjoint, hence: 



P(\ r n-n\ r ) 



exp ^- Leb 2 ^ |J ^.rl ^ exp (-£ Leb 2 (-Fb.r)) • 



To estimate this probability, let us introduce the triangle T r which is the 
greatest triangle included in Eq^\ \Li^^ l E k)r \ (see 3). We easily get: 

Leb 2 (T r ) = 2Rj8 + 0{\). 

For each k, E kr contains a triangle that is isometric to T r and disjoint from 
all others E k i r , hence we have: 

Leb 2 (\J E ki ^ > (£-1) Leb 2 (T r ) + Leb 2 (E , r ) , 

> £ Leb 2 (T r ) + Leb 2 (-E ,r) - Leb 2 (T r ) . 
Thus for all choice of £ > 2 and < &i < & 2 < • ■ ■ < kg < N r we have: 
exp(-Leb 2 (£ 2 (0,p r )))exp(-£Leb 2 (F , r )) < P (nf =1 F fc . )r ) < 
exp (- Leb 2 (£ ,r) + Leb 2 (T r )) exp (-^(Leb 2 (F 0>r ) + Leb 2 (T r ))) . 
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Fig 3. Left: definition of T r . Right: the set Eh x , T U ■ ■ ■ U Et e ,r in light grey and in dark 
grey the subset with area Leb2 (Eo,r) + (£ — 1) Leb2 (T r ). 



The number of such terms is ( N [), hence the sum Si of all those terms 
satisfies 

\S t \ < exp (Leb 2 (T r ) - Leb 2 (£ ,r)) exp (-£(Leb 2 (F , r ) + Leb 2 (T r ))) 

and the residual term S = X)£ 2 ( — 1)^ _1 <S^ is bounded from above by: 
\S\ < exp (Leb 2 (T r ) - Leb 2 (£ ,r)) x 
Nr fN\ 

J2[ I) exp (-^(Leb 2 (F ,r) + Leb 2 (T r ))) , 



£=2 



< exp (Leb 2 (T r ) — Leb 2 (Eq^)) x 

(1 + exp (-(Leb 2 (fb, r ) + Leb 2 (T r )))f r 



1 



N r exp (-(Leb 2 (F , r ) + Leb 2 (T r )))) , 



AT 2 

< -j- exp (Leb 2 (T r ) - Leb 2 (E , r )) x 

(5.7) exp (-2(Leb 2 (F , r ) + Leb 2 (T r ))) (1 + 0(1)) . 

Using the asymptotic expansions of Fo ir , Eq :T and T r and the unpper bound 
N r < C r ) we obtain: 

|S| < lexp^-^ KS + 0(l))cVx 

exp (-2R (2 - -k) r + 0(1)) (1 + 0(1)) 



< CCV exp ( -2Rr ( 1 + ( 1 - -k 



3 
I' 
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hence S = o(£r exp(— 2Rr)), which completes the proof of estimate 3.2. 

Remark 5.3. In the proof above, we could have taken only one tri- 
angle to obtain that the sum S is negligible with respect to the first term 
£s exp (— 2i?r), however the accuracy of the development would have been 
less interesting. Discs with bounded random radius R G [R*, R*] can also be 
treated this way, at a cost of a loss on the accuracy because of a non-optimal 
size of the fingers. 

Appendix 2: proof of Lemma 3.1. For sake of simplicity, we call 
x = ru. Let us denote by y (resp. y') a point in the intersection of (ru + L) 
with its tangent line emanating from O and situated on the left-hand side 
(resp. on the right-hand side) of the half-line (O + R+u). We define z (resp. 
z') as the orthogonal projection of y (resp. y') on {O + R+u) and a (resp. 
a') as the angle between (O + R+u) and (O + R+y) (resp. (O+R+y')). Then 
we have 

/'111/ — ^ 1 1 \ f \ \ — ^ 1 1 \ 

(5.8) $(ru + L) = a + a' = arctan ^ — — - — J + arctan y — — / ' 

Let us now describe W U (L): there exist two points w and w' (w being on 
the left-hand side of (O + M+u)) such that 

(5.9) W U (L) = dist(w, O + Ru) + dist(w / , O + Ru) 

where dist(-, O + Ru) is the Euclidean distance to the line (O + Ru). Com- 
paring (5.8) with (5.9), we observe that we only have to prove that 

( I \y — ^11^ 

(5.10) lim ra = lim r arctan I — — - — I = dist(w,0 + Ru) 



-oo r — ' + oo 



and 



f\W — z'\\\ 

(5.11) lim ra' = lim r arctan — — — — = dist^', O + Ru), 

the limits being uniform as required. Let us now concentrate on the first 
limit (the second can be proved in the same way): 

since ||y — z\\ < \\y — x\\ < D and \ \z — x\\ < \\y — x\\ < D, we have for 
every r > D 

I \v — z\ I | \y — z\ I D 

(5.12) tan a = < < -. 

\\z\\ r — L> r — D 
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Moreover, 

< dist(u>, O + Ru) - | \y - z\ \ = \ \y - w\ \ sm((3) < D sin(/3) 

where (3 is the angle between (O + M + u) and the line from y to w. Since y 
is a contact point of a support line of L and w is in L, this angle (3 must 
necessarily be lesser than a. Consequently, we get by a direct use of (5.12) 
that 

D 2 

(5.13) < dist(w,0 +lu) - \\y - z\\ < Dsin(a) < Dtan(a) < 



r-D 

Inserting this last estimate in the first equality of (5.12), we have 

) - 1)2 

r + D 



_ dist(u;, O + Ru) - \ /dist(^«,0 + Mu) , 

r arctan < ra < r arctan 



/dist(io,0 + Mu)\ 
V r-D J 



which provides the required convergence result (5.10) with a uniformity in 
u and in L. 
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